Cuprates as doped ^(1) spin liquids 
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We explore theoretically the notion that the underdoped cuprates may be viewed as doped (7(1) 
spin liquid Mott insulators. We pursue a conceptually clear version of this idea that naturally incor- 
porates several aspects of the phenomenology of the cuprates. We argue that the low doping region 
may be fruitfully discussed in terms of the universal physics associated with a chemical potential 
tuned Mott transition between a U{1) spin liquid insulator and a d-wave superconductor. A precise 
characterization of the deconfinement in the U{V) spin liquid is provided by the emergence of a con- 
served gauge flux. This extra conservation law should hold at least approximately in the underdoped 
materials. Experiments that could possibly detect this conserved gauge flux are proposed. 



I. INTRODUCTION 

Superconductivity in the cuprate materials occurs 
upon doping Mott insulators. The evolution of the phys- 
ical properties from the insulator to the superconduc- 
tor as a function of doping has been extensively stud- 
ied in the last several years. Many experiments have 
clearly established the relevance of the proximity to the 
Mott state at zero doping for understanding phenomena 
at finite doping where superconductivity appears at low 
temperature^. 

There is considerable theoretical debate however on 
the precise nature of the connection between the un- 
doped Mott state and the doped superconductor. Exper- 
imentally it is well established that the Mott insulating 
state at zero doping develops long ranged Neel antifer- 
romagnetism. Upon doping the magnetism disappears 
rather rapidly and soon thereafter is replaced by the su- 
perconducting ground state. The "normal" state above 
the superconducting transition temperature is metallic 
and is commonly referred to as the pseudogap state. A 
key observation is that the pseudogap state and the su- 
perconductor that descends from it both remember the 
Mottness of their parent undoped material but do not 
remember its long ranged Neel order. This observation 
is encoded in the old suggestions*^ that a useful way to 
think about the underdoped region is to first understand 
the nature of the possible non-magnetic Mott states at 
zero doping - obtained for instance by increasing frus- 
tration in the interaction between the spins. The be- 
haviour of the doped system may then possibly be fruit- 
fully viewed as the result of doping this non-magnetic 
Mott state. The idea is that the doping effectively frus- 
trates the Neel order so that the system is pushed across 
the transition where the Neel order is lost (see Fig. QJ. 

In this paper we pursue a particular non-magnetic 
Mott state that is connected to the Neel state by a 
second order transition. Furthermore doping this non- 
magnetic state leads to a d-wave superconductor with 
nodal quasiparticle excitations. For these and other rea- 
sons discussed below, the particular possibility we explore 
is theoretically very appealing. We show that this route 
to superconductivity leads to very unusual (and hence 
possibly unique) physical effects that could possibly be 



detected in experiments. 

The particular non-magnetic state we consider may be 
dubbed a 'deconfined' U{1) spin liquid on the two dimen- 
sional square lattice. The low energy theory of this state 
consists of nodal linear dispersing Dirac spinous that are 
coupled to a non-compact fluctuating U{1) gauge field. 
Such a spin liquid state has played a central role in a large 
number of prior theoretical papers in the field^'^'^'^. Re- 
cent developments^ have considerably clarified the con- 
ceptual basis of this state opening the way to the perspec- 
tive used in this paper. Indeed the present paper may be 
viewed as providing a conceptually clear version of vari- 
ous previous theoretical ideas that allows for predictions 
that are possibly testable. 

What precisely is meant by the notion that the un- 
derdoped cuprates may be viewed as doped versions of 
a non-magnetic Mott insulator? To understand this it is 
instructive to consider the phase diagram as a function 
of the chemical potential rather than the hole doping 
as shown in Fig. |21 Consider any non-magnetic Mott 
state that when doped leads to a d-wave superconduc- 
tor. As a function of chemical potential, there will then 
be a zero temperature phase transition where the holes 
first enter the system. For concretenes we will simply 
refer to this as the Mott transition - the corresponding 
phase boundary is marked in red in Fig |21 The associ- 
ated quantum critical fixed point will control the physics 
in a finite non-zero range of parameters. The various 
crossovers expected near such transitions are well-known 
and are shown in Fig. O Sufficiently close to this zero 
temperature critical point many aspects of the physics 
will be universal. The regime in which such universal be- 
havior is observed will be limited by 'cut-offs' determined 
by microscopic parameters. In particular we may expect 
that the cutoff scale is provided by an energy of a frac- 
tion of J (the exchange energy for the spins in the Mott 
insulator). We note that this corresponds to a reasonably 
high temperature scale. 

Now consider an underdoped cuprate material at fixed 
doping X. Upon increasing the temperature this will fol- 
low a path in FigO that is shown schematically. The 
properties of the system along this path may be usefully 
discussed in terms of the various crossover regimes in Fig. 
In particular it is clear that the 'normal' state above the 
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FIG. 1: Schematic zero temperature phase diagram showing 
the route between the antiferromagnetic Mott insulator and 
the d-wave superconductor. The vertical axis is labelled by a 
parameter g which may be taken as a measure of the frustra- 
tion in the interaction between the spins in the Mott insulator. 
AF represents the antiferomagnetically ordered state. SL is 
a spin liquid insulator that could potentially be reached by 
increasing the frustration. The path taken by the cuprate ma- 
terials as a function of doping x is shown in a thick dashed-dot 
line. The question marks represent regions where the physics 
is not clear at present. Doping the spin liquid naturally leads 
to the dSC state. The idea behind the spin liquid approach is 
to regard the superconducting system at non-zero x as result- 
ing from doping the spin liquid though this is not the path 
actually taken by the material. 
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FIG. 2: Same as in Fig. but as a function of chemical 
potential rather than hole doping. 



superconducting transition is to be understood directly 
as the finite temperature 'quantum critical' region asso- 
ciated with the Mott transition. Empirically this region 
corresponds to the pseudogap regime. Thus our assertion 
is that the pseudogap regime is controlled by the unstable 
zero temperature fixed point associated with the (Mott) 
transition to a Mott insulator. 

Once we adopt this point of view it is clear that the 



FIG. 3: Schematic phase diagram for a doping induced Mott 
transition between a spin liquid insulator and a d-wave super- 
conductor. The bold dot-dashed line is the path taken by a 
system at hole density x that has a superconducting ground 
state. The region marked FS represents the fluctuation regime 
of the superconducting transition. The region marked QC is 
the quantum critical region associated with the Mott critical 
point. This region may be identified with the high tempera- 
ture pseudogap phase in the experiments. 



Mott transition (which controls the pseudogap regime) 
need not be to the antiferromagnetic Mott state real- 
ized in the parent materials. Indeed the physical system 
takes a path through a complicated region which may in- 
volve charge ordering or stripe phases whose relation to 
other regions of the phase diagram is poor;y understood. 
Rather for the reasons mentioned in previous paragraphs 
it makes sense to explore Mott transitions to insulators 
that are non-magnetic. In the bulk of this paper we will 
focus on the U{1) spin liquid insulator and the Mott tran- 
sition from it to a superconductor induced by doping. 
Later we show how many aspects of the resulting theory 
may still survive provided only certain weaker assump- 
tions are satisfied (see Section rvij) . 



II. MOTT INSULATOR 

We begin with a discussion of the Mott insulator. The 
Neel state is the ground state of the nearest neighbour 
Heisenberg antiferromagnet. However as explained above 
we will be interested in possible non-magnetic ground 
states that are proximate to the Neel state. It is very 
important to realize that there are very few such known 
candidate non-magnetic ground states. Thus this theo- 
retical approach to the cuprate problem is fairly tightly 
constrained. One natural candidate is the dimcrized 
state described in Ref. ^ Studies of the doped dimcr- 
ized state have been pursued with some phenemenolog- 
ical successiS. Indeed a superconducting state with d- 
wave symmetry obtains upon doping. However such a 
superconductor also inherits the dimer order of the par- 
ent Mott state. In particular it has a full gap to spin 



excitations (at least at low doping) and breaks transla- 
tion symmetry. Empirically however there is strong evi- 
dence for the presence of such nodal quasiparticles, and 
for the presence of translation symmetry in most of the 
materials studied. We are thus naturally led to search 
for translation invariant non-magnetic Mott states that 
when doped will produce a d-wave superconductor with 
nodal quasiparticles. 

Mott insulators that preserve translation and other 
lattice symmetries are rather exotic beasts - the exci- 
tation spectra of all known theoretical examples are con- 
veniently described in terms of fractionalized spin-1/2 
'spinon' degrees of freedom. Further in certain such 'spin 
liquid' states, it is possible for the spinous to have gap- 
less nodal points with linear dispersion. Such spin liquids 
therefore provide ideal starting points to dope to produce 
translation invariant d-wave superconductors with nodal 
quasiparticles. 

An important physical property of such fractionalized 
spin liquid Mott insulators is the emergence of extra 
topological structure. Specifically spin liquid states pos- 
sess extra topological conservation laws not present in 
the microscopic models in which they arise. This con- 
servation law is conveniently interpreted as the flux of 
an emergent gauge field. Different classes of spin liquids 
may be defined depending on the nature of this conserved 
gauge flux. 

A simple example of a spin liquid with nodal fcrmionic 
spinous is provided by a state where there is a con- 
served Zi gauge flux. Indeed an attractive theory of 
the cuprates based on doping this state was developed 
and advocated in Ref>ii*i^. One feature of this theory 
is that in the superconducting state hcje vortices tend 
to have lower energy than than hc/2e vortices particu- 
larly at low doping. This is because in this route to su- 
perconductivity a hc/2e vortex necessarily involves the 
presence of Z2 gauge flux in its core. There is an en- 
ergy cost to having the Z2 gauge flux which dominates 
in the low doping limit and raises the energy of the hc/2e 
vortex. A crucial experimental teslii^ of such a theory is 
to directly detect the 'vison' excitation that carries this 
conserved Z2 flux, or indirectly to look for signatures of 
stable hc/e vortice a^^i^^ . Unfortunately, to date, all such 
experiments have been negativei^^i^. 

This forces us to reflne our search further and look for 
spin liquid states that do not possess Z2 gauge structure 
(i.e do not have vison excitations). An immediate candi- 
date for such a state is suggested by slave particle mean 
field theories of spin-1/2 models. Within a fermionic rep- 
resentation of the spin, a popular mean field state is the 
d-wave or staggered fiux state. This state is conveniently 
viewed as a d-wave paired state of spinous at 1/2 filing 
with the following Hamiltonian: 

H = - Y. Xrr'ifUr' + h.c) + ^rr'{frxfr'i - (t^i))- 

(2.1) 
Here /^q, a =],[ is a two-component fermionic spinon 
on the site r of a two dimensional square lattice. The 



hopping Xrr' is a real constant, and the pairing A^r' 
has d-wave structure: A^r' = Aq on horizontal bonds 
and Arr' = — Aq on vertical bonds. Diagonalization of 
the mean field Hamiltonian gives a low energy spectrum 
with 4 nodal points near which a linear dispersing Dirac 
spectrum appears. 

As pointed out in Refiii, this mean field state is equiv- 
alent under a unitary transformation to the 'staggered 
fiux' Hamiltonian: 



H 



sf 



= J2 {iiXrr'+Arr')Prfr'+h.c) (2.2) 



Here we take r to belong to one sublattice of the square 
lattice so that r' belongs to the opposite sublattice. This 
describes fermionic spin-1/2 spinous on the square lat- 
tice with complex hopping amplitudes such that there is 
a non-zero fiux that is staggered from plaquette to pla- 
quette. 

The crucial question is the fate of this mean field state 
upon including fluctuations. It is well-known that the im- 
portant fluctuations involve coupling to a compact U{1) 
gauge field. The resulting Hamiltonian is conveniently 
written in the 'staggered flux' gauge and has the struc- 
ture: 

H = Hf+Hg (2.3) 

Hf - 5I(*^-' + Ar.')e""-' Prfr' + h.c) (2.4) 

rr' 

Hg = -K^cosiV xa)+u ^ e^^, (2.5) 

P <rr'> 

Here a^r' G [0, 2tt) is to be regarded the spatial compo- 
nent of a U{1) gauge field, and Crr' is the corresponding 
electric field. Strictly speaking we must take the limit 
K —>■ 0,u -^ oo but relaxing this condition is not ex- 
pected to crucially change the physics. The Hamilto- 
nian must be supplemented with a Gauss law constraint 
V.e + flfr = 1 on every lattice site. 

Specializing to the low energy limit, this theory 
may formally be viewed as a theory of massless Dirac 
fermions coupled to a compact U{1) gauge field. Here as 
usual compactness means that point-like instantons or 
monopoles are allowed in the configurations of the gauge 
field in space-time. At each such monopole event the 
gauge fiux changes by an integer multiple of 27r. Recent 
work* has resolved a long-standing controversy on the 
low energy behavior of theories of this type, and shown 
that when the number TV of Dirac species is large, there 
is a stable 'deconflned' phase where the instanton fugac- 
ity renormalizes to zero at long distances and low ener- 
gies. (The case of direct physical interest corresponds to 
iV = 4) . A continuum fleld theory that describes this low 
energy flxed point is simply given by the action: 
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fer 



d2a;dT*„7^(a^ -1- M^)*a -I- -^{e^^xd^axf 

^^9 



(2.6) 



where ^^ , j = 1 , . . . . , 4 represent the four species of Dirac 
fermion appropriate for each of the two distinct nodes. 

A precise characterization of the deconfinemenljSii^ in 
this phase is obtained by noting that the absence of 
instantons at low energies imphes that the gauge flux 
b — dxCLy — dyUx is conserved. This corresponds to an ex- 
tra global topological U{1) symmetry at the low energy 
fixed point that is absent in the microscopic spin model. 

The analysis of Ref. |3 shows that stable U{1) spin 
liquids could exist in d = 2 at least for large enough N . 
Whether this stability extends down to the physically 
relevant case of SU(2) spins is not known at present. 
This corresponds to the case TV = 4. However the large- 
N analysis shows that as a matter of principle there is 
no reason to dismiss the possibility of such stable U{1) 
spin liquids for SU{2) spins. For most of the present 
paper we will simply assume that this is true, and explore 
its consequences. However detailed numerical calculation 
to settle this issue will certainly be useful and welcome. 
Towards the end of the paper (in Section IVip we will 
explore the opposite possibility - that the stability of the 
[/(I) spin liquids does not extend to SU{2) spins. We 
show how even in that case it may be possible to retain 
the key aspects of the theory provided certain weaker 
assumptions are satisfied. 

This d-wave paired U{1) spin liquid state is connected 
to a conventional Neel state at ordering vector (tt, tt) by a 
second order transition. This transition will be discussed 
elsewhere and is an example of a deconfined quantum 
critical poinljiSiiS. In particular there will be two diverg- 
ing length scales on the ordered side - one characterizing 
the crossover from critical to Neel ordered spin corre- 
lations, and a different longer one associated with the 
confinement of spinous. 

A schematic phase diagram depicting the Neel and spin 
liquid phases at zero doping is shown in Fig. ^ Now 
consider doping the Mott insulator. As in many previ- 
ous works, we assume that the doping moves the system 
across the phase boundary where the magnetism is lost. 
We may then fruitfully view the doped non-magnetic 
state as the result of doping the C/(l) spin liquid. 



III. DOPED [/(I) SPIN LIQUIDS 

What then happens to the U{1) spin liquid when it 
is doped? One answer to this question is suggested b 
the slave particle mean field theory developed in Ref. 
which also provides useful mathematical formalism. A 
d-wave superconductor with gapless nodal quasiparticles 
is obtained. To understand this it is first useful to ask 
about the nature of the electric charge carrying excita- 
tions in this doped U{1) spin liquid. Speaking loosely 
in the presence of fermionic spin-1/2 electrically neutral 
spinous, the charge of the doped holes may be expected 
to be associated with spin-0, charge-e bosons which carry 
gauge charge. This is roughly correct - more precise con- 
sideration shows that there are two species of such bosons 



which each carry gauge charge ±1. These will be denoted 
&i, 62. The presence of two species of bosons with gauge 
charges ±1 (but electric charge e and spin-0) is actually 
common in situations that involve doping U{1) fraction- 
alized phases. For instance, in three dimensional frac- 
tionalized boson insulators with a deconfined U{1) gauge 
field such as that considered in Ref. l20l two boson species 
appear quite naturally. Note that a composite of the two 
bosons is gauge invariant, has charge 2e, and is a spin 
singlet. Hence it may be identified with a Cooper pair. 
Thus the two bosons may each be identified as being one 
half of a Cooper pair. 

Now consider condensing each of the two species of 
bosons with the same amplitude < 61 >=< 62 >^ 0. 
Clearly this state breaks the physical electromagnetic 
U{1) symmetry, and hence is superconducting. Further 
once the bosons condense the spinous get endowed with 
electric charge - the structure of the spinon Hamiltonian 
imposes d-wave symmetry on this superconducting state. 
In particular the low energy nodal spinons simply become 
the nodal quasiparticles of the superconducting state. 
This condensate with < 5i >=< 62 >^ is actually 
energetically favored in slave particle mean field calcula- 
tions on the t — J model at low temperature in the low 
doping range. 

As discussed in the introduction it is instructive to 
consider the phase diagram not as a function of doping 
X but as a function of the chemical potential. The phase 
diagram then looks as shown in Fig|3| As the chemi- 
cal potential is increased, there is a Mott transition from 
the U{1) spin liquid insulator to the d-wave supercon- 
ductor. The corresponding (unstable) zero temperature 
fixed point controls many aspects of the physics in the 
underdoped side. In particular it determines the univer- 
sal aspects of the phsyics of the pseudogap regime and 
its eventual low temperature transition into the super- 
conducting state. 

What is the theory that describes this Mott transi- 
tion? To understand this we first note that in the insu- 
lating state, the two bosons 61^2 are gapped. Increasing 
chemical potential decreases this gap. The gap closes at 
the Mott transition beyond which superconductivity is 
achieved when both boson species condense (with equal 
amplitude). A continuum theory for the Mott transition 
that describes this condensation is readily written down. 
It takes the form 



S — Sfer + Sbi + Sb2 
Sbl = 



(3.1) 
bl[dr-tao-fi- ^^~ '"^ ) 61 (3.2) 



Sb2 = / b2 \ dr +iao- ^I- 



2mb 

(V + zd)'^ 
2mb 



62 (3.3) 



The change in the sign of the terms involving the gauge 
potential (ag , a) reflects the opposite gauge charges car- 
ried by 61 and 62 . The physical hole density is simply the 
total boson density given by 6J61 + 62^2, and thus cou- 



pies linearly to the chemical potential /i. The fermion 
part Sff^r is given in Eqn. 12.61 above. The transition 
occurs when ^ = 0. As discussed in the introduction a 
material at fixed doping x will take a path such as that 
shown in Fig|2|in this phase diagram. 

In writing down this continuum theory we have as- 
sumed that the instantons are irrelevant not just at the 
[/(I) spin liquid fixed point but also at the critical fixed 
point describing the Mott transition. This is quite rea- 
sonable. The presence of extra gapless matter fields is 
only expected to drive the instantons more irrelevant. 
With this assumption the conservation of gauge flux con- 
tinues to hold right at this critical fixed point. 

In Ref. la, a slave-particle representation with a 
fermionic charge-0 spin-1/2 spinon and two species of 
bosonic charge e spin-0 holons was introduced for the 
t-J model. Clearly this representation is ideally suited to 
discuss the zero doping U{1) spin liquid and the corre- 
sponding doped system. In particular the slave particles 
are 'deconfined' in the U{1) spin liquid and the associ- 
ated Mott critical point to the doped superconductor in 
the sense that instantons in the associated U{1) gauge 
field are irrelevant. As emphasized above this leads to an 
extra global topological symmetry associated with con- 
servation of gauge flux. This conserved gauge flux 
gives precise meaning to the notion of deconfine- 
ment and justifies the use of the slave particles as 
the useful degrees of freedom. Indeed the consider- 
ations of the previous paragraphs give sharp meaning to 
the theory of the underdoped cuprates developed on the 
basis of this slave particle representation. 



IV. MEANING OF GAUGE FLUX 

As the conserved gauge flux plays a crucial role in 
this theory it is interesting to ask what its meaning is 
and if its presence can directly be detected in the non- 
superconducting state just above Tc in underdoped sam- 
ples. In the following we describe some physical effects 
that should obtain if such a deconfined C/(l) gauge field 
indeed exists. The key is to exploit the unusual structure 
of the hc/2e vortex that obtains in the superconducting 
state. The unusual nature of the vortex stems from the 
fact that the superconducting state is obtained by con- 
densing separately the two charge-e bosons bi^2- Never- 
theless the superconducting state is smoothly connected 
to a regular BCS d-wave superconductor. Thus, despite 
the condensation of charge-e bosons, the superconductor 
must support hc/2e vortices. Such vortex solutions were 
described by Lee and WenSi. A simple construction in 
fact gives two apparently distinct hc/2e vortices - in one 
there is a full 27r vortex in 61 but not in 62 and in the 
other there is a 27r vortex in 62 but not in bi. Note that 
the Cooper pair operator ~ &1&2 winds by 2tt in either 
case. To understand these vortex solutions better and as 
they will prove important for what follows, we consider 
the following simple energy functional which captures the 



basic physics associated with the vortices: 



E = / d'x- 



V01 



A 



ve-i 



A 



(4.1) 
Here 6*1^2 are the phases of 61.2 respectively. The field a 
is the vector potential of the internal C/(l) gauge field, 
and A is that of the external electromagnetc field. The 
ellipses represent various other terms (such as the gauge 
and external field kinetic energies, etc) that have not been 
written out. Let ^1,2 wind by 27rTOi,2 on going on a large 
loop encircling a vortex (toi,2 are integers). Then we 
must have 



27rmi = <j)a + 4>A 
27rm2 = -0a -I- 0A 



(4.2) 
(4.3) 



Here (t>a.A are the fluxes of the gauge fields a, A enclosed 
by the loop. For a /ic/2e vortex we must have (pA = tt. 
This will happen if mi -t- 7n2 = 1 . The internal gauge flux 
<f>a = 7r(?ni — 7712) can then be any odd integral multiple 
of TT. In general there wil be some energy cost to having 
non-zero internal gauge flux so the lowest value of the 
gauge flux is energetically preferred. This still leaves the 
two possibilities (j>a = ±7r. These correspond precisely to 
the two vortices described above. Ref. |211 also pointed 
out that as the amplitude of only one of the two bosons 
will be suppressed in the core of either vortex, the region 
in the core looks like a condensate of 61 but not 62 or 
vice veras. Such condensates can easily be seen to have 
true staggered flux order with physical currents circulat- 
ing in a staggered manner in the core. This staggered 
current pattern corresponds to a broken symmetry phase 
that apparently rears its head in the vortex core. The 
symmetry of translation by one lattice spacing is broken 
by the staggering of the current pattern. However as the 
core is flnite the symmetry cannot truly be broken, and 
there will be a finite tunneling rate between the two kinds 
of cores. 

It is clear from the discussion above that the two kinds 
of cores and the associated staggerd patterns correspond 
to 4>a = ii". We immediately see that the tunneling event 
between the two kinds of cores must involve a change of 
2tt in the flux of the internal gauge field. Such a 27r flux 
changing event is precisely the instanton in the internal 
gauge field. Thus we have the remarkably simple picture 
of the instanton as just the tunneling event between the 
two kinds of staggered flux cores for hc/2e vortices in the 
superconducting state. 

Note that by assumption the superconducting state de- 
scends from a 'normal' state with deconfinement where 
the instantons are irrelevant at long scales. However in 
the superconducting state the internal gauge flux associ- 
ated with the vortex is confined to a flux tube of finite 
size. The finite size of the flux tube implies that in- 
stantons can no longer necessarily be ignored. Indeed 
instantons will also render unstable vortices that have 

(pa — ±37r, ±57r, so that there is in principle a unique 

/ic/2e vortex. 



These vortices will presumably start developing in- 
tegrity in the fluctuation region (see Fig|3Jl above the 
superconducting transition. For materials with small Tc 
{i.e close to the Mott critical point in Fig.), the instanton 
rate will be small (compared to electronic energy scales). 
Consequently even in the fluctuation region the vortex 
cores will have staggered flux order that fluctuates slowly 
at this rate. It is natural to identify this region with the 
one observed to have an enhanced Nernst effect in the 
experiments^. 

The considerations on vortices above enable us to pro- 
vide some meaning to the gauge flux once there is rea- 
sonably well-developed local superconducting order {i.e 
in the fluctuation region or below). In this region, a use- 
ful description is provided by a 'sigma-model' approach 
which focuses exclusively on the fluctuating classical or- 
der parameters. Indeed such a description was obtained 
in Ref. zj. Here we will take a slightly different perspec- 
tive that will bring out some key physical aspects. 

The fluctuations of the superconducting order parame- 
ter 6{62 *re clearly important in this region. As discussed 
above the vortices in this order parameter will start hav- 
ing integrity in this region. Since as we have seen above 
the vortices have slowly fluctuating staggered flux order 
in their cores it will be necessary to also include an order 
parameter for the staggered flux fluctuations. Clearly 
this is just simply given by n^ ~ |6ip — |62p- Note that 
both these order parameters (superconducting and stag- 
gered flux) have simple descriptions as bilinears in the 
underlying bosons. 

It is tempting to combine these order parameters into a 
single three-component vector n with n+ being the pair- 
ing order parameter and n^ the staggered flux one. In 
the regime where there is well-developed local order one 
might imagine that this vector has well-deflned magni- 
tude but there are fluctuations in its direction. A de- 
scription in terms of a fluctuating unit vector fleld which 
represents the local ordering direction may then possi- 
bly be appropriate. The energy functional for this unit 
vector fleld will presumably have easy plane anisotropy 
which favors the superconducting order. 

Despite the appeal of such a description it actually 
hides a very important piece of physics, and hence is not 
the full story. The point is that in the original description 
in terms of the boson flelds 61^2, there is an extra con- 
served quantity - namely the gauge flux - which should 
be included in the long wavelength physics. This extra 
conservation arises due to the assumed irrelevance of in- 
stantons at the Mott critical fixed point. (Of course at 
finite temperatures close to but not quite at the critical 
point the instanton fugacity will be small and hence strict 
conservation of gauge flux will not obtain - nevertheless it 
is necessary to include it to meaningfully discuss physics 
in the region impacted by the proximity to this critical 
point). 

What is this conserved flux in terms of the unit vector 
fleld n? The answer to this is well-known: it is the den- 
sity of defect configurations known as skyrmions. Thus 



in the absence of instantons, the skyrmion number in 
the 'sigma'-model description is conserved. It is impor- 
tant to appreciate that this is not a property of familiar 
quantum 0(3) sigma models in two dimensions but is 
rather closely related to the unconventional 0(3) models 
studied in Refs. 1241251129 . The latter are models in which 
hedgehog configurations of the unit vector fleld have been 
artificially suppressed by hand. In this situation as in the 
present problem a gauge theoretic description in terms of 
CP^ fields is superior as it directly brings out the extra 
conservation law. 

An expression for the conserved gauge flux in terms of 
the observable order parameter fields is readily written 
down. The gauge magnetic field 



B = 2Trpsk 



h.dxfi X dyh. 



We write 



fi = (sin cos a. sin sin a, cosij 



(4.4) 
(4.5) 

(4.6) 



so that n^ — cos 6 is the staggered flux order param- 
eter, and a is the phase of the superconducting order 
parametep^ Then the expression above for the gauge 
magnetic field is readily manipulated to 



i3 = i ( Vn^ X Va 



(4.7) 



It is easy to check that for a vortex with the structure 
described above the total gauge flux is ±n as expected. 



V. DETECTION OF GAUGE FLUX 

What are some signatures of the gauge flux in experi- 
ments? The asymptotic conservation of the gauge flux at 
the Mott transition flxed point potentially provides some 
possibilities for its detection. At non-zero temperatures 
in the non-superconducting regions, the flux conservation 
is only approximate (as the instanton fugacity is small 
but non-zero). Nevertheless at low enough temperature 
the conserved flux will propagate diffusively over a long 
range of length and time scales. Thus there should be an 
extra diffusive mode that is present at low temperatures 
in the non-superconducting state. It is however not clear 
how to design a probe that will couple to this diffusive 
mode at present. 

Alternately the vortex structure described above pro- 
vides a useful way to create and then detect the gauge 
flux in the non-superconducting normal state. We will 
first describe this by ignoring the instantons completely 
in the normal state. The effects of instantons will then 
be discussed. 

Consider first a large disc of cuprate material which 
is such that the doping level changes as a function of 
the radial distance from the center as shown in Fig. 0] 
The outermost annulus has the largest doping xi. The 




FIG. 4: Structure of the sample needed for the proposed ex- 
periment. The outer annulus (in dark blue) has the highest 
Tc. The inner annulus (in light blue) has a smaller Tc- The 
rest of the sample (in brown) has even smaller Tc- 



inner annulus has a lower doping level X2. The rest of 
the sample is at a doping level x^ < X2 < xi. The 
corresponding transition temperatures T'ci,2,3 will be such 
that Tc3 < Tc2 < Tel. We also imagine that the thickness 
ARo,ARi of the outer and inner annuli are both much 
smaller than the penetration depth for the physical vector 
potential A. The penetration depth of the internal gauge 
field a is expected to be small and we expect it will be 
smaller than Ai?o, Ai?^. We also imagine that the radius 
of this inner annulus Ri is a substantial fraction of the 
radius Ro of the outer annulus. 

Now consider the following set of operations on such a 
sample. 

(i) First cool in a magnetic field to a temperature Tin 
such that Tc2 < Tin < Tci- The outer ring will then 
go superconducting while the rest of the sample stays 
normal. In the presence of the field the outer ring will 
condense into a state in which there is a net vorticity on 
going around the ring. We will be interested in the case 
where this net vorticity is an odd multiple of the basic 
hc/2e vortex. If as assumed the physical penetration 
depth is much bigger than the thickness ARg then the 
physical magnetic fiux enclosed by the ring will not be 
quantized. 

(ii) Now consider turning off the external magetic field. 
The vortex present in the outer superconducting ring will 




FIG. 5: Variation of the doping level x with radial distance r 
in the sample. 



stay (manifested as a small circulating persistent current) 
and will give rise to a small magnetic field. As explained 
above if the vorticity is odd, then it must be associated 
with a flux of the internal gauge field that is ±7r. This 
internal gauge flux must essentially all be in the inner 
'normal' region of the sample with very small penetration 
into the outer superconducting ring. It will spread out 
essentially evenly over the full inner region. We have 
thus managed to create a configuration with a non-zero 
internal gauge flux in the non-superconducting state. 

(iii) How do we detect the presence of this internal 
gauge flux? For that imagine now cooling the sample fur- 
ther to a temperature T/i„ such that Tcs < ?/;„ < Tc2- 
Then the inner ring will also go superconducting. This 
is to be understood as the condensation of the two boson 
species 61,2. But this condensation occurs in the presence 
of some internal gauge flux. When the bosons 61,2 con- 
dense in the inner ring, they will do so in a manner that 
quantizes the internal gauge flux enclosed by this inner 
ring into an integer multiple of tt. If as assumed the inner 
radius is a substantial fraction of the outer radius then 
the net internal gauge flux will prefer the quantized val- 
ues ±7r rather than be zero (see below) . However config- 
urations of the inner ring that enclose quantized internal 
gauge fiux of ±7r also necessarily contain a physical vor- 
tex that is an odd multiple of hc/2e. With the thickness 
of the inner ring being smaller than the physical pen- 
etration depth, most of the physical magnetic flux will 
escape. There will still be a small residual physical flux 
due to the current in the inner ring associated with the 
induced vortex. This residual physical magnetic flux can 
then be detected. 

Note that the sign of the induced physical flux is inde- 
pendent of the sign of the initial magnetic field. Further- 
more the effect obtains only if the initial vorticity in the 
outer ring is odd. If on the other hand the initial vor- 
ticity is even the associated internal gauge flux is zero. 
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and there will be no induced physical flux when the inner 
ring goes superconducting. 

The preceding discussion ignores any effects of instan- 
tons. In contrast to a bulk vortex in the superconducting 
state the vortices in the set-up above have macroscopic 
cores. The internal gauge flux is therefore distributed 
over a region of macroscopic size. Consequently if in- 
stantons are irrelevant at long scales in the normal state, 
their rate may be expected to be small. At any non-zero 
temperature (as in the proposed experiment) there will 
be a non-zero instanton rate which will be small for small 
temperature. 

When such instantons are allowed then the internal 
gauge flux created in the sample after step (ii) will fluc- 
tuate between the values -I-tt and — tt. However so long 
as the time required to form the physical vortex in step 
(iii) is much longer than this instanton rate we expect 
that the effect will be seen. 



A. Energy estimates 

In this subsection we take a closer look at the ener- 
gies involved with various conceivable outcomes of the 
process described above. We assume that the cooling is 
slow enough that the system always stays in equillibrium 
so that the outcome of the experiment is determined by 
thermodynamic considerations. This will give us some 
indication of the conditions under which the experiment 
has its best chance of giving a positive result. We will 
use the simple energy functional in Eqn. I4.1l but will now 
need to keep track of the kinetic energy terms for the in- 
ternal gauge field. In principle we should also account 
for the existence of many weakly coupled layers in the 
three dimensional sample and the energy of the external 
magnetic field. To begin with we will ignore these latter 
effects and comment on them later. 

The energy per layer in the superconducting state is 
thus given by 



E 



Eh + En 



(5.1) 



h'4\{ 



Eh = / d^x— 



\79i-a-A 



Ea 



d xg ( V X a 



(5.3) 



Here K is the superfluid stiffness and is roughly of the or- 
der of ksTc. The parameter g is more difficult to estimate 
in any reliable manner. A crude estimate is provided by 
the following consideration. Right at the Mott critical 
point, we may hope to calculate the energy associated 
with various configurations of the gauge field in a 'RPA' 
approach which integrates out both the spinon and boson 
fields. In this calculation the gauge action is dominated 
by contributions from the spinous and is roughly given 
by 



'e//[ 



[a] r^ Ja d^qdujy'uj^ + q^\a^{q,uj)f (5.4) 



with J the exchange energy and a the lattice spacing 
(not to be confused with the vector potential). For the 
energy, this corresponds to an effective coupling g that is 
singular: g{q) ~ Ja/q for small q. At finite temperature 
T, this singular behaviour will be cutoff at a value qx ~ 
kBT/vp ^ kBT/Ja. We thus have the estimate 



5^ 



keT' 



(5.5) 



We caution however that this estimate is very crude. 

The energy associated with a supercurrent (the K 
term) competes with the energy associated with the in- 
ternal gauge field a. At the end of the experiment dis- 
cussed above, the system has one of two options. It can 
choose to quantize the gauge flux in the inner ring at 
TT while paying the cost of a spontaneous supercurrent 
through the inner annulus. The other option is that it 
prferes to push all the internal gauge flux out into the 
region between the two superconducting annuli thereby 
creating no supercurrent. If this were to occur the ex- 
periment will give a null result. It is therefore necessary 
to compare the energies of these two possible outcomes. 

The energy cost of the supercurrent is readily esti- 
mated and is 



Ej - {AR,)R, 



K 



Tc 



AR, 

~r7 



(5.6) 



For a sample with Tc ~ lOK, and ARi ~ O.lRi, we have 
Ej ^ IK. The energy cost of pushing the fleld out of 
the inner ring is roughly 



Ea--g 



1 



2„2 



J^a 



Rl - R} keTRlf 



(5.7) 



where /ttR^ is the area between the inner and outer rings. 
If we assume / ~ 0.1, then Ea > IK for samples with 
Ro <~ 1000a. Thus the energy cost per layer of the 
gauge field may be expected to dominate for samples of 
outer diameter about or smaller than a micron. However 
in view of the crudeness of the energy estimates, this 
should only be taken as very rough guidance. 

It is clearly also advantageous to have several layers 
so that the total energy gain in forming the supercurrent 
is multiplied by the number of layers and can exceed 
ksT. Finally we also mention that if the inner annulus 
thickness is much bigger than the physical penetration 
depth then all of the physical fiux associated with the 
trapped inner vortex will stay in the sample. However 
the energy cost of this physical fiux is easily seen to be 
enormous - it overwhelmes any effects due to the internal 
gauge field. It is for this reason that we advocate using 
an inner annulus thinner than the physical penetration 
depth to let most of the fiux escape. 



VI. ALTERNATE THEORETICAL POSSIBILIES 

In this Section, we wish to sketch some alternate the- 
oretical possibilities that weaken the assumptions made 



in previous Sections. In particular let us contemplate the 
situation that the stability of the U{1) spin liquid state 
found in the limit of a large number N of species fo Dirac 
spinons does not extend to the physically relevant case 
of A^ = 4. In other words, what if such a C/(l) spin liq- 
uid is never a stable phase for any model of SU{2) spins 
with short ranged interactions in two dimensions? If the 
instability is due to the relevance of a monopole opera- 
tor then it is natural to expect confinement to occur. A 
natural possible confined state that can result is simply 
the conventional Neel antiferromagnet^. 

We now show how even in this case it may still be 
possible to view the underdoped cuprates as doped U{1) 
spin liquids for the physics in an intermediate window 
of temperatures. This window will describe the pseudo- 
gap regime. Our arguments will rely crucially on lessons 
drawn from recent w orh^^^ on deconfined quantum crit- 
icality. Indeed Ref. Il8ll9l showed that unstable gapped 
U{1) spin liquids may emerge as good descriptions of the 
physics over a broad intermediate region of length and 
time scales close to certain quantum phase transitions. 
Here we will argue that in a similar vein the U{1) spin 
liquid state with gapless Dirac spinons considered in this 
paper may also control the physics in a broad interme- 
diate regime near certain quantum transitions at zero 
doping even if it is eventually unstable toward confine- 
ment. We will then discuss the implications for the doped 
system. 

We now consider the possibility that the spin liquid 
state in Figs. ^ 121 is a Z2 state of the spin system that 
has nodal Dirac spinons. We will argue that a large re- 
gion near the transition between this state and the con- 
ventional Neel state may be controlled by the (unstable) 
U{1) spin liquid fixed point. We recall that the Z2 state 
also has a gapped 'vison' excitation that carries the Z2 
gauge flux. Now consider a confinement transition out of 
this state that is obtained by condensing the vison. For 
gapped Zi spin liquids, such a transition to a confined 
valence bond solid state is possible. We refer the reader 
to Ref. [19! for a complete discussion. As discussed there, 
this transition provides an example of a deconfined quan- 
tum critical point. In particular the critical theory is that 
of a critical charged boson (interpreted as a spinon pair 
field) coupled to a fluctuating non-compact U(\) gauge 
field. 



To study the confinement transition out of the Z2 spin 
liquid with gapless nodal spinons, we follow the same 
general idea as in Ref. "lOl Such a state is conveniently 
described as a theory of Dirac spinons coupled to a com- 
pact C/(l) gauge field provided we also condense a singlet 
pair of spinons. Such a spinon pair has gauge charge-2 so 
that when it condenses, the [/(!) gauge theory enters a 
Higgs phase whose universal physics is described in terms 
of a deconfined Zi gauge theory. Schematically consider 



an action with the following general structure 

S = Sfer + Sfffi + Ssp (6-1) 

5/0 = / d^xdrn [(f>spip cr^ip + h.c) (6.2) 

d^xdT\{d^ - 2ia^)(j)sp\^ + .... (6.3) 



Here (psp represents the spinon pair field - this is taken 
to be invariant under all the global symmetries of the 
microscopic spin model. As mentioned above the Z2 spin 
liquid appears as the Higgs phase where (jjgp is condensed. 
Now consider a phase transition where 4)sp uncondenses, 
i.e gets gapped. To discuss the low energy physics of 
such a phase we simply drop terms involving (psp. The 
result is precisely the theory of Dirac spinons coupled 
to a compact U{1) gauge field. In this Section we are 
assuming that monopoles are relevant at the U{\) spin 
liquid fixed point for the physical case iV = 4 (TV is the 
number of Dirac species) . The resulting state is expected 
to be confined, presumably just the usual Neel state. 

It is however perfectly possible that monopole events 
though relevant at the U{\) spin liquid fixed point are 
irrelevant at the critical fixed point controlling the tran- 
sition to the Higgs phase {i.e the Z2 spin liquid). In- 
deed this is precisely what happens in the situations with 
gapped spinons considered in Ref. llMig . A particu- 
larly close analog is provided by the transition between 
a valence bond solid paramagnet and a spin gapped Z2 
paramagnet. Furthermore in the present gapless case, 
the critical fixed point has a greater number of gapless 
matter fields than in the U(l) spin liquid fixed point - in- 
creasing the number of gapless matter fields is expected 
to increase the scaling dimension of the monopoles. Thus 
it appears possible that monopoles are irrelevant at the 
critical fixed point though relevant at the f/(l) spin liq- 
uid fixed point. In the following we will assume that this 
is the case. This is a weaker assumption than the one 
made in previous sections. 

Within this weaker assumption, a number of inter- 
esting possibilities arise. First it strongly suggests the 
possibility of a direct second order transition between 
a Z2 spin liquid with Dirac spinons and a conventional 
collinear Neel antiferromagnet. The critical point is de- 
scribed by the theory in Eqns. 16.11 and where the com- 
pactness may be ignored at low energies. This is a de- 
confi ned qu antum critical point in the same sense as in 
Ref. Il8ll9i In particular the flux of the C/(l) gauge field 
is conserved at the low energy fixed point. The assumed 
renormalization group flows near this critical fixed point 
are shown in Fig. As in other situations with these 
flows, on approaching the transition from the Neel side, 
there will be two diverging length scales S, and S,conf with 
^conf diverging as a power of ^. Thus £,conf is potentially 
much bigger than ^ near the transition. The first length 
scale ^ is the scale below which the correlations appear to 
be that of the critical fixed point. At length scales bigger 
then ^ but smaller than ^conf, a description in terms of 
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U(l) SL 



FIG. 6: Possible renormalization group flow diagram that 
allows for a direct seconf order transition between a conven- 
tional coUinear Neel magnet, and a Z-i spin liquid with gap less 
Dirac spinous, g is the parameter used to tune through the 
transition. A is the fugacity of the leading allowed monopole 
operator. The (7(1) spin liquid with gapless Dirac spinous 
appears as an unstable fixed point that nevertheless controls 
the physics in a broad intermediate window of length scales 
in the ordered side of the phase transition. 



\ 



'conf 



Critical 



U(l) SL 



Neel 



FIG. 7: Crossovers as a function of length scale in the Neel 
state close to the transition assuming the validity of the RG 
flows of Fig. 1^ There is a broad intermediate regime where 
the physics is that of gapless Dirac spinous in a UiV) spin 
liquid state. 



the C/(l) spin liquid with Dirac spinons is appropriate. 
Finally it is only on length scales bigger than i^conf does 
the system crossover to the confined Neel state. Thus 
though unstable the t/(l) spin liquid would (within the 
assumptions above) still describe the physics in a broad 
intermediate region near the transition between the Neel 
state and a Z2 spin liquid with Dirac spinons. 

What may we then say about the doped system? It 
is now only necessary to assume that doping effectively 
pushes the system toward the critical point to the Z2 
spin liquid. Then it is natural to expect that the doped 
version of the U{1) spin liquid provides a description of 
the physics in a broad intermediate regime of length and 
energy scales. This intermediate regime may possibly 
describe the pseudogap region in moderately doped sam- 
ples. 

At present we do not know which (if cither) of the two 
possibilities for the fate of the U{1) spin liquid discussed 
above obtains. We hope that the considerations here will 
set the stage for future work on these questions. 



VII. CONCLUSION 

We have pursued the old idea that the underdoped 
metallic cuprates are fruitfully viewed as doped non- 
magnetic Mott insulators. We sharpened the meaning 
of this notion by first considering the phase diagram as 
a function of chemical potential rather than hole density. 
We suggested that the doped metal may be understood as 
being close to a chemical potential tuned quantum tran- 
sition from a spin liquid Mott insulator to the d-wave su- 
perconductor. From this perspective the actual detailed 
path followed by the real material as it evolves from the 
undoped antiferromagnetic Mott insulator to the doped 
metal is not that important in determining many aspects 
of the latter. Indeed this path in the real material ap- 
pears to be complicated and probably non-universal (in 
that it varies between different cuprates). What is uni- 
versal is the behavior in a wide window of intermediate 
temperatures (or energy scale) in the doped metal. Our 
suggestion is that this physics is captured by the univer- 
sal properties of the Mott critical point between the spin 
liquid insulator and the d-wave superconductor. 

The particular spin liquid state we considered has gap- 
less nodal spinons coupled to a fluctuating non-compact 
C/(l) gauge field. These spinons may be considered 'de- 
confined'. We emphasize that this term does not mean 
that the spinons behave as free fcrmions because they 
are still strongly coupled to the fluctuating C/(l) gauge 
fieldi. Rather the precise meaning to the notion of de- 
confinement is provided by the observation that there is 
an extra emergent global (topological) conservation law 
in this phase that is not present at the microscopic level. 
This conserved quantity is simply the gauge fiux. We 
discussed the meaning of the gauge flux, and suggested 
physical effects that could help detect its existence in ex- 
periments. 
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